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LONG-TERM  GOALS 

Wind-wave  interactions  are  of  fundamental  importance  as  they  determine  the  sea  surface  drag  and 
scalar  exchange  between  the  atmosphere  and  ocean.  This  is  particularly  important  at  high  winds  since 
air-sea  coupling  controls  tropical  cyclone  (hurricane)  formation  and  intensity.  There  is  large 
uncertainty  in  the  bulk  momentum  and  heat  surface  exchange  coefficients  (Cd,  Ck)  derived  from  field 
observations  and  laboratory  experiments  for  varying  wind  speed,  wave  age,  swell  amplitude  and 
direction,  and  in  the  presence  of  spray,  with  an  even  greater  debate  as  to  the  underlying  dynamical 
processes  that  couple  the  winds  and  waves  (Black  et  al.,  2007,  French  et  al.,  2007,  Powell  et  al.,  2003, 
Sullivan  et  al.,  2008).  In  this  context,  the  overarching  theme  of  this  project  is  to  identify  the  physical 
processes  that  couple  winds  in  the  marine  atmospheric  boundary  layer  and  the  underlying  surface 
gravity  wave  field  at  high  wind  speeds. 

OBJECTIVES 

This  project  seeks  to  reconcile  laboratory  and  field  measurements  of  wind-wave  interaction  and 
surface  drag  in  high  to  extreme  winds  using  turbulence  resolving  large-eddy  simulation  (LES).  The 
basic  science  question  we  address  is:  how  confidently  can  we  upscale  dynamical  processes  and 
measured  statistics  in  small-scale  laboratory  experiments  to  full-scale  high  wind  ocean  conditions? 

APPROACH 

We  build  substantially  on  our  accumulated  expertise  in  sea  surface  processes  and  air-sea  interaction. 


1 


Report  Documentation  Page 


Form  Approved 
OMB  No.  0704-0188 


Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 
VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 


1.  REPORT  DATE 

30  SEP  2014 


2.  REPORT  TYPE 


3.  DATES  COVERED 

00-00-2014  to  00-00-2014 


4.  TITLE  AND  SUBTITLE 

Turbulence  Simulation  of  Laboratory  Wind-Wave  Interaction  in  High 
Winds  and  Upscaling  to  Ocean  Conditions 

6.  AUTHOR(S) 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

National  Center  for  Atmospheric  Research, Boulder, Co 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 


5a.  CONTRACT  NUMBER 


5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

5d.  PROJECT  NUMBER 


5e.  TASK  NUMBER 


5f.  WORK  UNIT  NUMBER 

8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


10.  SPONSOR/MONITOR'S  ACRONYM(S) 

11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 


12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 


13.  SUPPLEMENTARY  NOTES 


14.  ABSTRACT 


15.  SUBJECT  TERMS 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 
ABSTRACT 

18.  NUMBER 
OF  PAGES 

19a.  NAME  OF 
RESPONSIBLE  PERSON 

a.  REPORT 

unclassified 

b.  ABSTRACT 

unclassified 

c.  THIS  PAGE 

unclassified 

Same  as 
Report  (SAR) 

7 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


For  this  project,  LES  of  the  airflow  over  steep  and  breaking  waves  is  carried  out  over  a  range  of  wind 
speeds  and  quasi-steady  wave  train  conditions  matching  selected  laboratory  wind-wave  experiments 
(Banner,  1990).  More  recently  we  extended  our  study  to  include  transient  steep  water  wave  packets 
that  we  computed  using  an  accurate  fully  nonlinear  boundary  element  wave  model  (Banner  et  al., 
2014).  Simulation  data  are  analyzed  extensively  to  generate  flow  visualization,  bulk  mean  flow  and 
turbulence  statistics  and  surface  drag,  and  are  compared  with  the  observed  experimental  results.  This 
detailed  intercomparison  of  simulation  and  experimental  results  is  used  to  identify  the  fundamental 
dynamical  processes  in  wind-wave  interaction,  viz.,  flow  separation  and  reattachment,  impacts  of 
surface  roughness  and  swell  steepness  as  a  function  of  non-breaking  and  breaking  conditions.  The 
second  stage  of  the  project  will  attempt  to  upscale  the  laboratory  dynamical  processes  to  full-scale 
oceanic  conditions  under  high  wind  and  hurricane  forcing  using  LES  of  the  full  marine  ABL. 

The  turbulence  simulation  method  uses  a  novel  approach  that  allows  coupling  turbulent  winds  to  an 
underlying  three-dimensional  time-dependent  wave  field.  In  the  simulations,  the  surface  wave  fields, 
which  contain  the  essential  features  of  breaking,  are  externally  prescribed,  i.e.,  space-time 
measurements  of  the  wave  fields  are  imposed  as  a  lower  boundary  condition  in  the  computations.  The 
LES  code  is  highly  parallelized  and  all  the  computations  are  performed  on  parallel  super-computers 
using  meshes  of  0(512  )  using  512  or  more  computational  cores. 

WORK  COMPLETED 

Our  effort  during  FY13  concentrated  on  refining  how  closely  the  LES  modeling  framework  developed 
by  Sullivan  can  replicate  the  findings  reported  in  the  laboratory  investigation  of  Banner  (1990),  which 
investigated  the  marginal  effect  of  breaking  on  the  airflow  over  a  quasi-steady  train  of  steep  wind 
waves  generated  by  a  constant  frequency  wave  paddle  and  a  constant  wind.  We  explored  the  sensitivity 
of  our  initial  LES  results  to  the  fluid  surface  speed  distribution  along  these  profiles,  wave  profile  shape 
and  background  roughness  length.  These  findings  were  documented  in  detail  in  our  FY13  report  and 
showed  close  qualitative  and  quantitative  correspondence  between  the  observed  and  modeled  physics. 

During  FY 14  our  efforts  focused  on  analyzing  and  refining  our  large-eddy  simulations  of  turbulent 
flow  over  incipient  and  active  breakers.  Details  of  the  computational  methodology  are  documented  in 
our  FY  13  report  and  in  a  new  publication  Sullivan  et  al.  (2014).  One  of  the  topics  we  wished  to  pursue 
further  in  FY  14  is  the  LES  solution  sensitivity  to  the  grid  mesh  resolution,  i.e.,  are  the  LES  solutions 
we  generate  for  flow  over  breakers  grid  independent.  In  previous  work,  we  found  that  some  LES 
statistics,  in  particular  high-order  moments,  are  sensitive  to  the  grid  mesh  resolution  for  boundary- 
layer  flows  with  flat  lower  boundaries  (e.g.,  Sullivan  and  Patton,  201 1).  Thus,  we  repeated  our  FY  13 
calculations  using  (512“x  128)  grid  points  compared  to  the  previous  computations  that  used  a  mesh  of 
(2562  x  128);  the  domain  size  (XL,  YL,  ZL)=( 5,  5,  1)  is  held  constant.  We  refer  to  these  LES  as  fine  and 
coarse  grid  solutions,  respectively.  The  new  fine  grid  computations  use  4  times  as  many  grid  points, 
but  require  8  times  more  computational  resources  to  achieve  the  same  number  of  large  eddy  turnover 
times  (~  50)  owing  to  the  reduction  in  the  timestep  on  the  fine  grid.  Roughly,  the  cost  of  each  fine 
mesh  calculation  is  375,000  core-hours  running  on  512  processors  on  a  Cray  XE6  at  the  Engineering 
Research  and  Development  Center  (ERDC).  We  find  the  new  fine  mesh  runs  are  similar  in  character  to 
the  previous  coarse  mesh  runs,  but  are  pleasingly  in  better  agreement  with  the  observed  wind  profiles 
(see  the  discussion  in  RESULTS).  Our  analysis  of  the  LES  solutions  also  includes  computation  of 
numerous  statistics,  extensive  flow  visualization,  and  animation  of  the  solution  fields.  We  also 
commenced  an  extension  of  the  work  to  make  the  waves  more  representative  of  natural  waves,  which 
typically  occur  as  unsteadily  evolving  groups.  We  have  commenced  simulations  in  which  we  replaced 
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the  fixed  steep  profile  uniform  wave  train  by  a  train  of  unsteady  steep  wave  packets  generated  by  our 
fully  nonlinear  wave  code  (see  Banner  al.  (2014)  for  details).  In  related  work,  an  archival  manuscript, 
now  in  press,  documents  our  LES  solution  algorithm  for  turbulent  over  a  spectrum  of  moving  waves  in 
a  stratified  marine  boundary  layer  (Sullivan  et  al.,  2014). 

RESULTS 

Sample  results  from  our  fine  mesh  simulations  over  incipient  and  active  breakers  are  provided  in 
Figs.  1-4.  The  visualization  in  Fig.  1  shows  contours  of  spanwise  vorticity  coy  =  du/dz  -dw/dx  for 
turbulent  over  an  incipient  breaker.  Here  we  choose  to  visualize  vorticity,  since  it  is  Galilean  invariant, 
i.e.,  it  is  independent  of  the  wave  phase  speed  c,  and  thus  is  one  of  the  useful  metrics  for  gauging  the 
presence  of  flow  separation  over  a  moving  boundary.  The  visualization  shows  an  intense  concentration 
of  vorticity  on  the  forward  face  of  the  wave  that  is  compacted  into  a  shallow  layer  at  the  wave  surface. 
Animations  further  show  that  this  vorticity  layer  intermittently  detaches  from  the  wave  surface  near  the 
wave  crest  and  then  often  oscillates  vertically  above  the  wave  trough  as  it  is  advected  downstream.  The 
vorticity  layer  serves  to  isolate  slow  moving  fluid  in  the  wave  trough  from  the  faster  moving  overlying 
turbulent  flow.  The  lower  panel  of  Figure  1  shows  that  the  vorticity  layer  near  the  surface  and  also  in 
the  separated  flow  region  contains  strong  vortical  cores;  these  vortices  are  identified  using  the  well- 
tested  and  documented  technique  described  by  Zhou  et  al.  (1999)  and  Adrian  (2007).  In  this  image 
notice  how  the  vortical  cores  leave  the  wave  surface  in  the  separated  regions  and  propagate  into  the 
outer  turbulent  flow.  Overall,  these  results  support  the  premise  that  turbulent  flow  separation  occurs 
over  incipient  and  active  breakers  in  agreement  with  the  observation  by  Veron  et  al.  (2007). 

In  our  FES  the  net  form  (or  pressure)  drag  D  is  the  spatial  integral  of  the  pressure-waveslope 
correlation;  in  the  FES  transformed  coordinates  (x,  y,  z)  =>  (C,  //,  0 

D=ffj£xd£dr)  (1) 

where  J  is  the  Jacobian  of  the  grid  transformation  between  physical  and  computational  space  and  QJJ 
is  essentially  the  waveslope.  Figure  2  shows  the  spatial  distribution  of  the  pressure-waveslope 
correlation  for  active  and  incipient  breakers.  Key  features  to  notice  are:  pCx/J  is  spatially  intermittent 
and  varies  considerably  with  the  spanwise  y-coordinatc  despite  the  two-dimensional  character  of  the 
underlying  imposed  waves.  Closer  inspection  of  Figure  2  also  shows  that  the  largest  values  of 
pressure-waveslope  correlation  frequently  occur  on  the  forward  face  of  the  wave,  i.e.,  where  positive 
pressure  is  well  correlated  with  positive  waveslope.  This  result  can  be  understood  by  noting  that  the 
vorticity  layer  in  Figure  1  separates  near  the  wave  crest  but  often  stagnates  (re-attaches)  on  the  forward 
face  of  the  downstream  wave  generating  high  local  positive  pressure.  There  are  also  significant 
contributions  to  the  total  drag  from  regions  just  forward  of  the  wave  crest  where  negative  pressure  is 
well  correlated  with  negative  waveslope.  The  local  value  of  the  pressure-waveslope  correlation  (the  red 
regions  in  Fig.  2)  often  exceed  the  mean  value  by  more  than  a  factor  of  5. 

Average  wind  profiles  from  the  computations  are  compared  with  the  observations  by  Banner  (1990)  in 
Fig.  3.  We  find  (very)  good  agreement  between  the  observations  and  the  fine  mesh  FES  solutions.  In 
this  comparison,  the  FES  winds  are  averaged  over  £  and  ;/  at  constant  values  of  the  wave-following 
coordinate  C  whereas  the  observations  are  collected  at  a  fixed  height  above  the  mean  water  level.  In 
the  FES  the  constant  C  surfaces  become  nearly  level  surfaces  at  the  vertical  location  where  the 
observations  are  collected.  Sullivan  et  al.  (2014)  shows  that  averaging  in  wave-following  coordinates 
is  advantageous  as  it  allows  computation  of  the  winds,  and  also  the  fluxes,  down  to  the  water  surface, 
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i.e.,  below  the  wave  crests.  We  speculate  that  the  improved  simulation-observation  agreement  is  a 
consequence  of  better  resolution  of  the  separation  and  re-attachment  points  (shown  in  Fig.  2),  and 
reduced  reliance  on  the  subgrid-scale  model.  At  constant  £V/ 1  the  leftward  shift  of  the  wind  profile  with 
active  breaking  compared  to  the  incipient  case  implies  a  greatly  increased  surface  roughness  z0- 


Models  of  the  wave-induced  boundary  layer  typically  adopt  a  constant  flux  assumption  near  the  water 
surface  and  then  seek  to  partition  the  momentum  flux  into  turbulent  and  wave-induced  pieces  as 
function  of  height  in  the  constant-flux  layer.  These  ideas  can  be  evaluated  using  our  LES  results  for 
flow  over  active  and  incipient  breakers.  The  steady  ensemble-average  streamwise  momentum  balance 
equation  in  our  LES  assuming  periodic  sidewalls  reads  (Sullivan  et  al.,  2014). 


\j)='M1j^+u(w'z')+Tfc- 


(2) 


where  the  large-scale  driving  pressure  gradient  Px  is  balanced  by  the  vertical  divergence  of  three 
terms,  viz.,  in  left-to-right  order,  pressure,  resolved  turbulence  and  subgrid-scale  (or  viscous)  terms.  In 
the  above  expression,  W  is  the  contravariant  flux  velocity  in  the  vertical  direction  normal  to  a  £  —  rj 


surface,  zt  is  the  grid  speed  of  each  point  in  the  computational  mesh,  and  <  .  >  denotes  a  horizontal 
average.  In  our  problem  the  lower  boundary  condition  is  W  =  Zt  at  C=  0  so  that  there  is  no  flow  across 
the  wavy  surface.  Fig.  4  shows  the  vertical  distribution  of  the  streamwise  momentum  flux  components 
contributing  to  the  right-hand-side  of  (2)  for  active  and  incipient  breakers.  First,  the  total  flux  varies 
almost  linearly  over  the  height  of  the  computational  domain  and  accurately  satisfies  global  momentum 
balance.  However,  when  the  fluxes  on  the  right-hand- side  of  (2)  are  plotted  on  a  logarithmic  vertical 
axis  as  in  Fig.  4  we  find  a  shallow  constant  flux  layer  near  the  water  surface  extending  up  to  about 
C/X  ~  0.06.  It  is  interesting  to  notice  that  this  is  approximately  the  maximum  crest  height  as  shown  in 
the  bottom  panel  of  Fig.  2.  In  our  computations  of  a  full  marine  boundary  layer  with  height  z,  ~  500  m, 
we  find  the  constant  flux  layer  is  deeper  extending  up  to  more  than  0. 1  C/zi.  In  the  present  calculations, 
the  pressure-waveslope  correlation,  turbulent  flux,  and  subgrid-scale  terms  all  vary  smoothly  with 
height  above  the  wave  surface,  and  it  is  readily  apparent  that  in  the  case  with  active  breaking  the 
pressure  drag  dominates  the  flux  at  the  water  surface  carrying  almost  75%  of  the  total  wind  stress. 


or 
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Fig.  1  Snapshot  of  instantaneous  spanwise  vorticity  coy  contours  in  an  x-z  plane  (left  panel)  from  LES  over 
incipient  breakers  using  a  fine  mesh  of  (5122x  128)  gridpoints.  The  vorticity  is  normalized  by  X/u*.  Notice 
a  vigorous  (positive-signed)  vorticity  layer  remains  attached  to  the  forward  face  of  the  wave  but  separates 
from  the  wavy  surface  beyond  the  wave  crest  and  propagates  into  the  overlying  turbulent  flow.  The  right 
panel  shows  that  this  vorticity  layer  contains  strong  vortices.  Here  the  vortical  cores  are  identified  by  the 
complex  eigenvalues  of  the  velocity  gradient  matrix  duj/dxj,  i.e.,  the  so-called  “Xcl  method”  of  Zhou  et  al. 
(1999).  The  vorticity  layer  leaving  the  crest  of  the  wave  is  suggestive  of  flow  separation  as  observed  by 
Veron  et  al.  (2007)  in  a  wind-wave  tank. 
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Fig.  2  Instantaneous  contours  of  pressure-waveslope  correlation  (or  form  drag)  pdh/dx  in  an  x-y  plane  for 
active  (left  panel)  and  incipient  (right  panel)  breaking.  For  reference,  the  two-dimensional  wave  height 
h(x,  t)  is  shown  in  the  bottom  of  each  panel.  Pressure  is  normalized  by  pu?~  (p  is  the  air  density),  and  the 
winds  and  waves  both  propagate  from  left  to  right.  Notice  a  large  fraction  of  the  total  form  drag  occurs  on 
the  windward  face  of  the  waveform,  i.e., where  positive  pressure  is  well  correlated  with  positive  waveslope. 
The  flow  separates  near  an  upstream  wave  crest  and  stagnates  (re-attaches)  on  the  windward  face  of  a 
downstream  wave  generating  high  positive  pressure.  The  form  drag  is  (74,  54)%  of  the  total  drag  for 
(active,  incipient)  breaking,  respectively.  The  red  contours  are  regions  where  the  local  pressure  drag  is 
more  than  5  times  the  mean  value  obtained  by  averaging  over  the  entire  domain. 


Fig.  3  Vertical  profile  of  average  streamwise  wind  <JJ>  normalized  by  surface  friction  velocity  u*  for 
incipient  (blue  lines)  and  active  breakers  (red  lines)  from  LES.  Solid  and  dash-dot  lines  use  grid  meshes  of 
(5122x128)  and  (2562  x  1 28)  gridpoints,  respectively.  A  reference  rough  surface  log-linear  wind  profile  for 
no  resolved  surface  waves  is  also  shown  by  a  solid  black  line.  The  vertical  offset  of  the  wind  profile  at 
constant  <U>/u *  reflects  the  increased  roughness  induced  by  breaking  waves.  The  spatial  averaging  in  the 
LES  is  carried  out  in  wave-following  coordinates  at  constant  f  values.  Solid  bullets  are  measurements  of 
mean  wind  speed  obtained  at  fixed  heights  measured  from  the  mean  water  surface  from  Banner  (1990). 
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U  momentum  flux  /  u2 

Fig.  4  Vertical  profiles  of  the  terms  in  the  ensemble  average  streamwise  momentum  budget  computed  in  the 
LES  wave-following  coordinates  (see  Sullivan  et  ah,  2014).  Flow  over  incipient  and  active  breakers  is 
indicated  by  dashed  and  solid  lines,  respectively.  The  contributions  to  the  streamwise  momentum  are: 
pressure-waveslope  correlation  (form  drag)  pCx/J  (red  curves),  vertical  turbulent  momentum  flux  u{W—zt ) 
(green  curves),  the  subgrid-scale  flux  (blue  curves)  and  their  sum,  the  total  flux  (black  curves).  Fluxes  are 
normalized  by  the  total  stress  u2  and  results  used  a  fine  grid  of  (5122  x  128)  points. 

IMPACT/APPLICATIONS 

This  study  is  aimed  at  advancing  fundamental  understanding  of  air-sea  interactions,  which  play  a 
critical  role  in  forecasting  marine  weather,  ocean  waves  and  upper  ocean  dynamics.  Central  to  air-sea 
interaction  mechanics  is  the  coupling  between  turbulence  in  the  surface  layers  of  the  marine  boundary 
layers  and  the  connecting  surface  gravity  wave  field,  which  determines  the  sea  surface  drag  and  the 
scalar  exchange  between  the  atmosphere  and  ocean.  The  major  impact  of  this  effort  will  be  a 
refinement  of  present  knowledge  of  the  influence  of  wave  breaking  and  complex  sea  surface 
topography  in  air-sea  interactions.  This  includes  a  reduction  in  the  present  large  uncertainty  in 
measured  values  of  the  momentum  and  scalar  exchange  coefficients  for  wind  speeds  greater  than  20 
ms  1  by  clarifying  their  dependence  on  numerous  physical  factors,  such  as  wind  speed,  wave  age,  wind- 
wave  direction,  wave  spectral  content,  flow  separation  and  surface  roughness. 
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